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Abstract 

We have measured the branching fraction, longitudinal polarization fraction f L , and the CP 
asymmetry coefficients A and S in B° — > p + p~ decays. These results are obtained from a 
253 fb -1 data sample containing 275 million BB pairs collected by the Belle detector run- 
ning at the KEKB e+e~ collider. We obtain B = [24.4 ± 2.2 (stat) (syst) ] x 1(T 6 , f L = 

0.951 to ol (stat) ±o:8ii ( svst )> A = °- 00 ±0.30 (stat) (syst), and S = 0.09 ±0.42 (stat) ± 
0.08 (syst). These values are used to determine the CKM phase angle 4> 2 via an isospin analysis; 
the central value and la error are (87 ± 17)°, and 59° < (p 2 < H5° at 90% CL. 

PACS numbers: 13.25.Hw, 11.30.Er, 12.15.Hh 
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The decay B° — > p + p~ has received much attention in the literature because it allows 
one to determine the CKM phase angle 4> 2 |l| with relatively little theoretical uncertainty. 
This is because the penguin contribution to B° — > p + p~ is small, as indicated by the small 
branching fraction for B° — > p°p° Q . The angle 2 is determined by measuring the decay 
time (At) distribution of B° and B decays; the difference between the distributions is fit to 
the function ^lcos(Am At) + iSsin(AmAt), where Am is the mass difference between the 
two B°-B° mass eigenstates, and the CP asymmetry coefficients A and S depend on 2 . To 
determine 4> 2 requires knowledge of the polarization of the p mesons, as different polarization 
amplitudes can give different values of A and S for the same 2 . For a negligible penguin 
contribution, the sign of S for the CP-odd transversity amplitude A ± will be opposite to 
that for the CP-even amplitudes A and An. In this paper we present a measurement of 
the branching fraction, longitudinal polarization fraction (f L ), and coefficients A and S for 
B°^p + p~ decays. 

The data sample consists of 253 fb _1 recorded by the Belle experiment running at the 
KEKB energy-asymmetric e + e~ collider The Belle detector Q consists of a silicon ver- 
tex detector (SVD), a 50-layer central drift chamber (CDC), an array of aerogel threshold 
Cherenkov counters (ACC), time-of-flight scintillation counters (TOF), and an electromag- 
netic calorimeter comprised of CsI(Tl) crystals. These detectors are located within a su- 
perconducting solenoid coil providing a 1.5 T magnetic field. An iron flux- return located 
outside the coil is instrumented to identify muons and detect mesons. Two inner detec- 
tor configurations were used: a 2.0 cm beampipe and a 3-layer silicon vertex detector were 
used for the first sample of 152 million BB pairs, and a 1.5 cm beampipe, a 4-layer silicon 
detector, and a small-cell inner drift chamber were used for the remaining 123 million BB 
pairs |5j. 

Candidate B° — > p + p~, p^ — > n^n decays are selected by requiring two oppositely- 
charged tracks satisfying p T > 0.10 GeV/c, dr < 2.0 mm, and \dz\ < 4.0 cm, where p T is 
the momentum transverse to the beam axis, and dr and dz are the radial and longitudinal 
distances, respectively, between the track and the beam crossing point. The tracks are fit 
to a common vertex, and the resulting \ 2 is required to be satisfactory. We require that 
the tracks be identified as pions based on information from the TOF and ACC systems, 
and from dE/dx measurement in the CDC. This information is combined into a likelihood 
for a track to be a pion (LA or a kaon (L K ). To identify pions, we require that the ratio 
L K /(L K + LA be less than 0.40. The efficiency of this requirement is approximately 89%, 
and the kaon misidentification rate is about 10%. Tracks are rejected if they satisfy an 
electron identification criterion based on information from the ECL system. 

The ir^ candidates are subsequently combined with 7r° candidates. The latter are recon- 
structed from pairs of photons having an invariant mass in the range 0.1178-0.1502 GeV, 
which corresponds to ±3<r in 7r° mass resolution. The energy of each photon in the labora- 
tory frame must be > 50 (90) MeV in the ECL barrel (endcap) region, which subtends the 
angular range with respect to the beamline of 32°-129° (17°-32° and 129°-150°). In order 
to reduce combinatorial background, n° candidates are required to have p > 0.35 GeV/c 
in the e + e~ center-of-mass (CM) frame. To identify p ± — > 7r ± 7r° decays, we require that 
m w ± n o be in the range 0.62-0.92 GeV/c 2 ; this window corresponds to ±2<r in the ^ir mass 
distribution. To further reduce combinatorial background, we require that each p candidate 
satisfy —0.80 < cos 6 < 0.98, where 9 is the angle between the direction of the tt° and the 
negative of the B° momentum in the p ± rest frame. 

To select B° — > p + p~ decays, we calculate the quantities M bc = \/El cam — p\ and 
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AE = E B — E hcam , where E B and p B are the reconstructed energy and momentum of 
the B candidate, and E heam is the beam energy, all evaluated in the CM frame. The AE 
distribution has a long tail on the lower side due to incomplete containment of the electro- 
magnetic shower in the ECL. We define a signal region 5.27 GeV/c 2 < M bc < 5.29 GeV/c 2 
and —0.12 GeV < AE < 0.08 GeV. The fraction of events having multiple candidates in 
this region is 9.5%. Many of these arise from fake 7r°'s combining with good tracks, and 
thus we choose the best candidate based on the mass difference |^ 77 — "\o|- From MC 
simulation we find that this criterion correctly identifies the B° — > p + p~ decay 87% (95%) 
of the time for longitudinal (transverse) polarization. Signal decays that have one or more 
daughters incorrectly identified but nonetheless pass all selection requirements are referred 
to as "self-cross-feed" (SCF) background. Most such decays have the 7r° daughter swapped 
with a 7T° originating from the rest of the event; this preserves the vertex position and thus 
the At value. 

We identify whether a B° or B° evolved and decayed to p + p~ by tagging the b flavor of the 
non-signal (opposite-side) B decay in the event. This is done using a tagging algorithm [f| 
that categorizes charged leptons, kaons, and A's found in the event. The algorithm returns 
two parameters: q, which equals +1 (—1) when the opposite-side B is likely to be a B° {B°); 
and r, which represents the quality of the tag as determined from MC simulation and varies 
from r = for no flavor discrimination to r — 1 for unambiguous flavor assignment. 

The dominant background is e + e~ — > qq [q = u, d, s, c) continuum production. We dis- 
criminate against this using event topology: e + e~ — > qq processes in the CM frame tend 
to be jet-like, while e + e~ — > BB events tend to be spherical. To quantify this sphericity, 
we calculate a set of 16 modified Fox- Wolfram moments and combine them into a Fisher 
discriminant J?}. We calculate a probability density function (PDF) based on this discrim- 
inant and multiply it by a PDF for cos 6^, where 9 B is the polar angle in the CM frame 
between the B direction and the positron beam direction. BB events are produced with 
a 1 — cos 2 8 B distribution while continuum events are produced uniformly in cos8 B . The 
PDFs for signal and continuum background are obtained using MC simulation and the 
data sideband 5.21 GeV/c 2 < M bc < 5.26 GeV/c 2 , respectively. We use the products of 
the PDFs to calculate a signal likelihood C s and continuum likelihood £ qq - and require that 
TZ = CJ (C s + C qq ) be above a threshold. This threshold is determined by optimizing a 

figure-of-merit N s / yj N s + N B , where N s (N B ) is the number of signal (background) events 
estimated to be in the signal region from MC simulation (extrapolating from an M bc side- 
band). For this calculation we assume a B° — > p + p~ branching fraction of 25 x 10~ 6 . As 
the tagging parameter r also discriminates against qq background, we divide the data into 
six r intervals (denoted by £ = 1 — 6) and determine the TZ threshold separately for each. 
The TZ requirement removes about 97% of continuum events while retaining about 62% of 
B° — > p + p~ signal. The efficiency of all selection criteria as determined from MC simulation 
is (2.19 ± 0.02)%. This value corresponds to f L = 1; the efficiency for our measured value 
of f L (which is slightly less than one - see below) is about 4% higher, and we include this 
difference as an additional systematic error in the branching fraction. 

We determine the signal yield in two steps: we first fit the M hc -AE distribution to obtain 
the fraction of B° — ► p + p~ + B° —>■ p ± 7r T 7r non-resonant decays (as this fit cannot distinguish 
between the two modes); we then fit the m^^o distribution to determine the non-resonant 
fraction and hence the p + p~ signal yield. 

The first fit is an unbinned maximum likelihood (ML) fit to the two-dimensional M bc -AE 
distribution in the wide range 5.21 GeV/c 2 < M bc < 5.29 GeV/c 2 and -0.20 GeV < AE < 
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0.30 GeV. The fit includes components for the signal and several backgrounds: continuum 
events, b — > c decays, and b — > u decays. The PDFs for signal and b —* u background are 
modeled by smoothed two-dimensional histograms obtained from large MC samples. The 
PDF for b-^c background is the product of a threshold ("ARGUS" [8|) function for M bc 
and a quadratic polynomial for AE, also obtained from MC simulation. The PDF for 
continuum background is taken to be an ARGUS function for M bc and a linear function for 
AE; the slope of the linear function depends on the tag quality (r) bin £, and all seven shape 
parameters are floated in the fit. The signal PDF is adjusted to account for small differences 
observed between data and MC for a high-statistics sample of B + — ► D°p + , D° — > K + 7t~tt 
decays, which, like the signal mode, contain two neutral pions. 

The b — ► u background is dominated by £?(°<+) —> pn, B — > a^, and B — ► a x p decays. 
As their contributions are small, their normalization is fixed to the yield obtained from MC 
simulation. For B + —^afn and B—>a 1 p, the branching fractions are unknown; we therefore 
estimate them to be 3 x 10 -6 and 2 x 10~ 6 , respectively, and vary these values by 50% and 
100%, respectively, to obtain the systematic error due to these estimates. There are thus nine 
free parameters in the fit, and the resulting fraction of B° —>p + p^ + non-resonant decays in 
the signal region is (3.0 ± 0.4)%. Figure ^ shows the final event sample and projections of 
the fit result. The \ 2 °f the projections divided by the number of bins is 0.97 for M bc and 
1.02 for AE. 

To distinguish B° — > p + p~ decays from non-resonant B° — > p ± 7r =F 7r and B° — > 7r + 7r°7r~7r° 
decays, we do an unbinned ML fit to the 7r ± 7r° mass distribution. For this fit we require 
5.27 GeV/c 2 < M bc < 5.29 GeV/c 2 and -0.12 GeV < AE < 0.08 GeV and fit the m^ 
distribution in the range 0.30-1.80 GeV/c 2 . Only one p candidate is required to satisfy 
0.62 GeV/c 2 < "nt^i^o < 0.92 GeV/c 2 ; the mass of the other p candidate is then fit. We 
include additional PDFs for non-resonant B — > pirn and B — > rrrr components; however, 
the fit result for the latter is negligibly small (<C 1%) and consistent with zero, and thus 
we set this fraction to zero. The PDFs for continuum and b — > c backgrounds are grouped 
together and taken from the data sideband 5.22 GeV/c 2 < M bc < 5.26 GeV/c 2 ; we use MC 
simulation to check that the shapes of these backgrounds and their ratio in the sideband 
region are very close to those in the signal region. We impose the constraint that the fraction 
of signal + non-resonant events in the m ±^ range 0.62-0.92 GeV/c 2 equals that which we 
obtained from the M hc -AE fit; there is then only one free parameter in the fit (f p7T7T ). The 
result of the fit is f pnT = (8.9 ± 6.2)%, and thus N pp = 142 ± 13 events. Figure |1 shows the 
data (two entries per event) and the projection of the fit result. The \ 2 °f the projection 
divided by the number of bins is 1.35. 

The branching fraction is evaluated as 

5(5°^V) = ; p NpP N , (1) 

where N pp is the number of B° — > p + p~ candidates (142 ± 13), -/V — is the number of 
BB pairs produced [(274.8 ± 2.3) x 10 6 ], e is the geometric acceptance and efficiency of 
selection criteria obtained from MC simulation [(2.19 ± 0.02)%], and e PlD is a correction 
factor for the pion identification requirement to account for small differences between data 
and MC (0.969 ± 0.012). Inserting these values into Eq. (0) gives a branching fraction of 
(24.4 ± 2.2) x 10~ 6 , where the error given is statistical. 

There are several sources of systematic uncertainty; their effect is evaluated by varying 
the relevant parameter (s) by ±lcr and noting the resulting change in the branching fraction. 
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FIG. 1: Left: projection in M hc for events satisfying —0.10 GeV < AE < 0.06 GeV. Right: projection 
in AE for events satisfying 5.273 GeV/c 2 < M hc < 5.290 GeV/c 2 . The top plots correspond to all 
r values, and the bottom plots correspond to 0.75 < r < 1.0. The curves show projections of the 
ML fit result: the dashed curve represents p + p~ + pirir events, the dotted curve represents qq 
background, the dot-dashed curve represents b — > c background, the small solid curve represents 
b— >u background, and the large solid curve shows the overall result. 



The different sources and the corresponding changes are: track reconstruction efficiency 
(1.2% per track); tt° reconstruction efficiency (4% per 7r°); the n ± identification efficiency 
e PID in Eq. (JTJ) (1.4%); MC statistics used to calculate the acceptance e in Eq. (JTJ) (1.0%); 
the dependence of the acceptance upon f L (+0.0%, —4.1%); the continuum suppression 
requirement (13%); and the fraction of signal + non-resonant decays used as a constraint 
in the rn n±w0 fit. This last error itself has several components: the calibration factors used 
to correct the M hc -AE distribution of MC B° — > p + p~ events (to better match the data); 
the M hc -AE shapes used for b — > c background; the fraction of b —>■ u background, which 
is obtained from MC simulation and includes several estimated branching fractions; and 
the AE range used in the fit. These components give an overall variation in f pp + f p7T7T of 
(+2.9%, —8.6%). We subsequently redo the rn^^o fit, varying f pp +f p7Tn by this amount; the 
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FIG. 2: The rn^+^a plus Tn^—^o distribution (two entries per event) for events satisfying 
5.27 GeV/c 2 <M bc <5.29 GeV/c 2 and -0.12 GeV <A£< 0.08 GeV. Superimposed are projections 
of the ML fit result. The dashed curve represents p + p~ signal; the dot-dashed curve represents 
non-resonant pirn decays; the dotted curve represents qq + (b—>c) + (b—>u) backgrounds; and the 
solid curve shows the overall result. 



resulting change in N pp is (+2.0%, —6.4%). Combining this in quadrature with the other 
systematic errors gives an overall systematic error of (+16%, —17%). Thus the final result 
for the branching fraction is 



B(B°^p + p') = [24.4 ±2.2 (stat)^(syst) 



x 10" 



(2) 



To determine the polarization of B° — > p + p~ decays, we do an unbinned ML fit to the 
helicity angle distribution F(cos9 1 , cos# 2 ), where 9 l (9 2 ) is the angle between the direction 
of the 7T° from the p + (p~) and the negative of the B° momentum in the p + (p~) rest frame. 
For a longitudinal polarization fraction f L , this distribution is 



f l cos 2 6 1 cos 2 6 2 + -(1 - f L ) sin 2 9 1 sin 2 e 2 



(3) 



In the fit, this PDF is multiplied by a two-dimensional acceptance function for (cos 6^, cos# 2 ) 
determined from MC simulation. The acceptance is modeled as the product ^(cos^)- 
A (cos 6> 2 ), where A is a polynomial function. 

We fit events satisfying 5.27 GeV/c 2 < M bc < 5.29 GeV/c 2 , -0.12 GeV <AE< 0.08 GeV, 
and 0.62 GeV/c 2 <m 7r±7r0 <0.92 GeV/c 2 . The likelihood function includes PDFs for signal, 
pixix non-resonant decays, and continuum, b— >c, and b— >u backgrounds. The PDF for non- 
resonant decays is taken to be constant, and the PDF for b^u background is taken from MC 
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FIG. 3: The cos6 1 plus 0036*2 distribution (two entries per event) for events satisfying 
5.27 GeV/c 2 < M hc < 5.29 GeV/c 2 , -0.12 GeV < AE < 0.08 GeV, and 0.62 GeV/c 2 < m n±n0 < 
0.92 GeV/c 2 . Superimposed are projections of the ML fit result. The dashed curve represents 
p + p~ signal; the dotted curve represents non-resonant pmr decays and qq + (b — » c) + (b — » u) 
backgrounds, and the solid curve shows the overall result. 

simulation. The PDFs for continuum and b— >c backgrounds are combined and determined 
from the data sideband 5.21 GeV/c 2 < M bc < 5.26 GeV/c 2 , -0.12 GeV < AE < 0.12 GeV; 
we use MC simulation to check that the shapes of these backgrounds and their ratio in 
the sideband region are very close those in the signal region. The fraction of signal + non- 
resonant decays is taken from the previous M hc -AE fit; the component f p7rn is taken from 
the previous ui ± a fit. The fraction of b — > u background is small and is taken from MC 
simulation. Since fr q q +h _> c \ = 1 — f pp — f p7Tn — f^ u , there is only one free parameter in the 
fit (f L ). The result of the fit is f L = 0.951 laolL where the error given is statistical. FigureEl 
shows the data and projections of the fit result. The x 2 °f the fit projections divided by the 
number of bins is 0.83 for cos 6^ and 1.05 for cos #2. 

For this measurement there are eight sources of systematic error; their effect is evalu- 
ated by varying the relevant parameter(s) by ±ler and noting the resulting change in f L . 
The different sources and the corresponding changes are: the signal + non-resonant fraction 
(+0.013,-0.012); the non-resonant fraction (+0.021,-0.020); the pion identification effi- 
ciency, which has a small effect upon the acceptance (+0.000,-0.004); misreconstructed 
B° — > p + p~ decays (+0.005,-0.000); the continuum suppression requirement (±0.013); 
possible interference with an L = tt^tt system produced in non-resonant p ir^n de- 
cays (±0.005); a very small bias in the fitting procedure measured from a large toy MC 
sample (+0.000,-0.005); and uncertainty in the continuum +(6 — > c) background shape 
(+0.004, —0.014). This last uncertainty is evaluated by taking the background shape from 
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alternative M bc and AE sidebands. Combining all systematic errors in quadrature gives 
an overall systematic error of (+0.029,-0.031). Thus the final result for the fraction of 
longitudinally polarized decays is 

h = 0-951 ± °;°i (stat) Iffi (syst) . (4) 

To fit for the CP asymmetry coefficients A and S, we divide the data into q = +1 
and q = — 1 tagged decays and do an unbinned ML fit to the respective At distributions. 
At KEKB, the T(4S) is produced with a Lorentz boost of /?7 = 0.425 nearly along the 
electron beamline (z); since the B° and B° mesons are approximately at rest in the e + e~ 
CM frame, At is determined from the displacement in z between the fcp and / tag decay 
vertices: At ps (z C p — z t!lg ) / [3^c. 

The likelihood function for event i is 

a = f$n^) PP + /&^(At)scF + f { xn^t)^ + (5) 

/ 6 ( l^(At) 6 _ + fti u V(At)^ u + fSv(At) qq , 

where the event weights /W are functions of M bc and AS and are normalized to the fractions 
of events obtained from the previous M hc -AE and m^±^ fits. The PDFs V(At) for b— >c 
and b— >u backgrounds are determined from MC simulation, and the PDF for continuum qq 
background is determined from an M hc sideband. We include an additional PDF for SCF 
background in which a 7r =l= daughter is swapped with a track from the rest of the event; this 
function is determined from MC simulation and is found to be exponential with an effective 
lifetime of about 1 ps. The weighting function / SCF is also determined from MC simulation; 
its normalization is 5.7% of all p + p~ candidates. 
The PDF V pp {At) is given by 

r+ao e H A *'l/ r s o r 

/ —r U-gAu; <(i) +g(l-2u; <(i) ) x (6) 

[^lcos(AmAt / ) + 5sin(AmAt / )]|i?(At (i) , At') dAt' , 

where R(At, At') is a resolution function determined from data, uj e is the mistag probability 
for the £th bin of the tagging parameter r, and Au e is a possible difference in uj £ between 
q = +1 and q — — 1 tags. The PDF V p7Tn is taken to be exponential with r = r B and is 
smeared by the same resolution function R. We determine A and S by maximizing the 
log-likelihood J2i log£j, where i runs over the 656 candidate events satisfying 5.27 GeV/c 2 < 
M bc <5.29 GeV/c 2 , -0.12 GeV<A£<0.08 GeV, and 0.62 GeV/c 2 <m^ <0.92 GeV/c 2 . 
The results of the fit are A = 0.00 ± 0.30 and S = 0.09 ± 0.42, where the errors given 
are statistical. These values are consistent with the no-CP- violation case ^4 = 5 = 0, and 
the errors are consistent with expectations based on MC studies. Figure |U shows the data 
and projections of the fit result for 0.0 <r <0.5 and 0.5 <r < 1.0 subsamples; the x 2 °f the 
projections divided by the number of bins is 0.98 and 0.97, respectively, which is satisfactory. 

For this measurement there are twelve sources of systematic error; their effect is evaluated 
by varying the relevant parameter (s) by ±lcr and noting the resulting change in A and S. 
The different sources and the corresponding changes are listed in Table |l] The error due to 
wrong-tag fractions is evaluated by varying u e and Au e values. The error due to b — >• c and 
continuum At distributions is evaluated by varying the ~40 parameters used to describe 
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FIG. 4: The At distribution for events satisfying 5.27 GeV/c 2 < M bc < 5.29 GeV/c 2 , -0.12 GeV< 
AS < 0.08 GeV, and 0.62 GeV/c 2 < m ± a < 0.92 GeV/c 2 . Superimposed are projections of the 
ML fit result, (a) q = +\ tags; (b) q = —l tags; (c) raw CP asymmetry for 0.0<r<0.5; (d) raw 
CP asymmetry for 0.5 <r< 1.0. For (a) and (b), the hatched distribution represents signal events 
and the dashed curve represents background. 



these distributions. The error due to the component fractions is evaluated by varying the 
fractions obtained from the M hc -AE fit. The error due to self-cross- feed (SCF) events is 
evaluated by varying both the SCF fraction (by 50%) and the effective lifetime used for SCF 
events (by 2a). The effect of a CP asymmetry in b— >c and qq backgrounds is evaluated by 
adding such asymmetries to the & — ^ c and qq At distributions; the sizes of these asymmetries 
are taken to be those obtained from fitting an M bc sideband in the data (these values are 
small and consistent with zero). The error due to a possible fitting bias is evaluated using 
a full MC simulation, and the error due to tag-side interference |9 ! ] is evaluated using a toy 
MC simulation. The error due to the transversity amplitudes A ± and A« (which may have 
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different values of A and S) is evaluated in two steps: first setting f L equal to its central 
value and varying the values of (A, S) for A ± and Au over the range —1 to +1; second, using 
common values of (A, S) for all transversity amplitudes (except for an extra minus sign for 
Aj_) and varying f L . Combining all systematic errors in quadrature gives overall systematic 
errors of (+0.10, —0.09) for A and ±0.08 for S. Thus the final result is 

A = 0.00 ± 0.30 (stat) to.H (syst) (7) 
S = 0.09 ± 0.42 (stat) ± 0.08 (syst) . (8) 



TABLE I: Systematic errors for the CP asymmetry coefficients A and S. 



Type 


8 A (xlO- 2 ) 


8S (xl0~ 2 ) 




+<x 


—a 


+<T 


—a 


Wrong tag fractions 


0.5 


0.6 


0.8 


0.8 


Parameters Am, t b0 


0.1 


0.1 


0.9 


0.9 


Resolution function 


1.3 


1.3 


1.3 


1.3 


Background At distributions 


1.6 


1.5 


2.3 


2.5 


Component fractions 


2.1 


2.6 


5.0 


4.5 


pTnr non-resonant fraction 


0.0 


0.0 


0.6 


0.6 


Self-cross-feed (SCF) fraction 


0.0 


0.0 


0.2 


0.2 


Background asymmetry 


0.0 


2.0 


0.0 


4.3 


Possible fitting bias 


0.0 


0.2 


0.6 


0.0 


Vertexing 


4.1 


2.8 


0.4 


0.4 


Tag-side interference 


7.2 


5.5 


0.3 


0.0 


fj_, /|| (transversely- 










polarized components) 


5.2 


4.5 


6.0 


4.7 


Total 


+10.2 


-8.6 


+8.4 


-8.4 



We use these values along with the measured branching fraction for B° —>■ p + p~ and 
previously-measured branching fractions for B + — > p + p° jiol ] and B° — > p°p° to constrain 
the angle 2 . We use the isospin relations of Ref. |ll[ (originally applied to the B — > nix 
system), neglecting a possible 1 = 1 contribution to the B°^p + p~ amplitude 0|. We first 
fit the measurements to obtain a minimum x 2 ; this is denoted Xmm- We then scan values 



of 02 from 0°-180° and for each value calculate the difference Ax = X 



We 



insert A% 2 into the cumulative distribution function for the x 2 distribution for one degree of 

11 cf) 2 values 
(87+17)°; 



freedom to obtain the confidence level (CL) for the <p 2 value. The curve 1— CL for all cf) 2 values 



is plotted in Fig. From this curve we obtain a central value and la error (f> 2 
the 90% CL interval around the central value is 59°<0 2 <115°. 

In summary, we have measured the branching fraction, longitudinal polarization fraction 
f L , and CP asymmetry coefficients A and S for B° — ► p + p~ decays. Our results are B = 



24.4 ±2.2 



h3.8 
-4.1 



x in-6 f — n QP,1 +0 033 +0.029 a 
x iu i JL — U.yOl _ .039 -0.031) 



0.00 ±0.30 and S 



0.09 ±0.42 ± 



0.08, where the first error listed is statistical and the second error listed is systematic. These 
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FIG. 5: Results of fitting the branching fractions B(B° -> p + p~), B{B + -» B(B°^p°p°), 
and the p + p _ asymmetry parameters .A and 5 for the angle (j) 2 , using the isospin relations of 
Ref. 11 1 . The vertical axis is one minus the confidence level (see text). The horizontal line at 
1 — CL = 0.317 corresponds to a 68.3% CL interval for </> 2 , and the horizontal line at 1 — CL = 0.10 
corresponds to a 90% CL interval. 



results are consistent with previously-published measurements |L3J. The result for f L is 
consistent with the theoretical expectation based on QCD factorization 14], and the results 
for A and S are consistent with the case of no CP violation. From an isospin analysis we 
obtain the constraint 59°<0 2 <115° at 90% CL. 
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